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Abstract 

The equation of state of saturated nuclear matter is derived using two different derivative- 
coupling Lagrangians. We show that both descriptions are equivalent and can be obtained 
from the a — lo model through an appropriate rescaling of the coupling constants. We 
introduce generalized forms of this rescaling to study the correlations amongst observables 
in infinite nuclear matter, in particular, the compressibility and the effective nucleon mass. 

PACS number(s): 21.65.+f, 21.30.+y, 97.60.Jd 
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1 Introduction 



The equation of state (EOS) of dense matter is known, experimentally, at just one density value 
Po = 0.15/fm 3 (symmetric nuclear matter) for a binding energy per particle E/A(po) = —16 
MeV. The compression modulus is more loosely determined ranging from as low as k = 100 
MeV to the highest estimate k = 344 MeV [|. The last piece of experimental information is 
the effective mass which, from the energy dependence of the proton-nucleus optical potential, 
can be constrained to lie in the range M* /M=0.Q — 0.9. 

Theoretically, the EOS is obtained using as starting point either non-relativistic interactions 
of the Skyrme type 0] , or relativistic mean- field models (RMF) of baryons and mesons 0, [| |], 
0. As an example of the non-relativistic calculations, the Sftm^-interaction [RJ, used successfully 
to reproduce ground-state nuclear properties, gives M*/M = 0.78 and a compressibility k = 217 
MeV. As for the RMF calculations, the original linear o — uj model |4| gives for the same 
observables, a smaller effective mass M* /M = 0.556 and a larger compressibility k = 540 MeV. 
The relativistic treatment of the EOS becomes necessary at extreme conditions of density and 
temperature such as those found in heavy-ion reactions at energies above 400 MeV/nucleon 
and type-II supernovas. In what follows we will concentrate on it. 

The linear a — uj (Walecka) model satisfactorily explains many properties of nuclear matter 
and finite nuclei with two free parameters. The resulting compressibility at saturation density, 
however, exceeds the experimental bound. A way out of this difficulty is to introduce non-linear 
scalar self-couplings |9|, |10[ . The resulting non-linear a — uj model reproduces well ground state 
nuclear properties and is renormalizable (though it shows some instabilities at high densities 
for low values of the compression modulus, k < 200 MeV). The number of free parameters in 
this case is four. 

One alternative approach which renders a satisfactory compression modulus without in- 
creasing the number of parameters is that advanced in Refs. O, [T2|, [TB| . It employs the same 
degrees of freedom and the same number of independent couplings present in the Walecka 
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model. The difference is that it introduces a non-renormalizable derivative coupling to the 
baryon field later adjusted to reproduce the experimental conditions at saturation. The results 
for the compression modulus and effective mass compare well with the Skyrme-type calcula- 
tions. The spin-orbit splitting in finite-nuclei, however, turns out to be smaller than required 
by the data. We shall refer to this description as Model I. 

A qualitative different model, though similar in spirit, obtains if instead of modifying the 
covariant derivative only the baryonic kinetic energy term is redefined. As in Model I, a suitable 
rescaling of the baryonic field leads to a Lagrangian describing baryons with effective mass M*. 
Such a model was originally suggested as a possibility by Zimanyi and Moskowski in Ref. 



(appendix) and implemented for a particular case by Delfino et al. [II]. The experimental 
nuclear matter saturation density and binding energy per nucleon are used to fit the two free 
parameters; the compressibility comes out rather small but within the bounds (156 MeV) and 
the spin-orbit splitting doubles the value obtained in Model I (though still short of the data). 
To this model we shall refer hereafter as Model II. 

The aim of this paper is to provide a unified discussion of models I and II. We will show 
that both are to be understood as generalized Walecka models where the meson couplings (g a 
and g^) become effective coupling constants. In Model I only the scalar coupling is affected; in 
Model II both scalar and vector couplings are modified. The discussion will lead us naturally to 
consider a generalization of the models above effected by the introduction of a third parameter 
a which rescales the effective coupling constants. We will explore the results for M*, k, and 
the mean-field vector (V) and scalar (S) potentials, within both models for different values of 
the parameter a. 

Our basic motivation is to sum up in a phenomenological meson-baryon model, with as few 
parameters as possible, the ability to reproduce the four nuclear-matter observables and to give 
a good description of ground-state properties of finite nuclei without the presence of instabilities 
at high densities. It is not a substitute to the more fundamental relativistic Brueckner-Hartree- 



Fock (RBHF) method [15, 0| [17| but a simplifying alternative to explore the nuclear matter 
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equation of state. 

It is important to stress that the same idea of density-dependent coupling constants underlies 
two recent approaches. The first one, known as relativistic density-dependent Hartree-Fock Jl8| . 
describes finite nuclei and nuclear matter saturation properties using coupling constants that 
are fitted, at each density value, to the RBHF self-energy terms. The good agreement obtained 
for the ground state properties of spherical nuclei lends support to a description involving 
these density-dependent coupling constants. In the second approach [T9|] , a similar density 
dependence of the parameters is achieved through projection of the meson-nucleon vertices into 
positive-energy space weighted with a parameter adjusted to reproduce saturation properties. 
Here also, the generalization of the model through the inclusion of adjustable parameters at the 
vertices permits to study the density dependence in an analytical manner. Both approaches 
have natural partners in the relativistic density- dependent Hartree-Fock of Li and Zhuo 
and in the Gmuca's model Our treatment of models I and II belongs to this category of 
descriptions involving density-dependent coupling constants. 

The outline of the paper is as follows: in the next section we present the two models, establish 
their form and relationship and give expressions for the generalized EOS and compressibility. 
In section 3 we present our results and discuss the correlation between the calculated nuclear 
matter observables. Finally, in the last section we present our conclusions. 



2 The Models 

We begin by introducing the non-linear Lagrangian density of Model I, 

C NL = $ {fl{ 9ui ) - m»M} i, - -F^F^ + ~m>^ + ~ (d.ad^a - m 2 y) , (1) 

where p{gu)) = l^Fl^g^) = ^{id^ — g u uj^) is the usual covariant derivative and the degrees of 
freedom are the baryon field ip, the scalar-meson field a, and the vector-meson field uA The 
real function m*(a) is unknown except for the fact that it must go to unity as the density 
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goes to zero and has to vanish at large densities where the effective mass approaches zero 
asymptotically. Indeed, the Dirac equation obtained from the Lagrangian density ([!]) gives 

■»» = -jjr. p) 

where M and M* are the bare and effective baryonic mass, respectively. 

We show next that performing a transformation on the spinor field Cnl can be derived from 
a Lagrangian density with derivative scalar coupling (DSC). The proposed Lagrangian density 



13] which generates C^l is given by 

C = ^ {K(<7)]-% w ) - M) V> - \F^F, U + l -mlu^ + \ fraPa - m> 2 ) . (3) 

Introducing the rescaled baryonic field ip — > [m*(<T)]^ip, we obtain from Eq.(|3]) the rescaled 
Lagrangian density 

C R = £ NL + $f, (4) 
where 3 is an imaginary contribution given by 

3=^7W0^1n(m>)). (5) 

This term does not carry any physical content and can be transformed away by substituting the 
baryonic kinetic energy term in Eq.fll]) for the symmetric derivative ^{t/jj^d^i/j — (d tJ "ip)'y ll if}}. 
The imaginary contribution S cancels after the field scaling. Once performed the field rescaling 
is equivalent to the replacement, 

{[mV)]- 1 ^) - M} - {p(g u ) - m»M} . (6) 

The Lagrangian densities given by Eqs.(^) and @ are, therefore, completely equivalent; they 
describe the same physics whether we deal with infinite nuclear matter or finite nuclei. The 
equations of motion obtained from Eqs.(0) and @ give rise to the same hadronic dynamics. 
The modified kinetic energy of Eq.(§) is not arbitrary; it describes the motion of a baryon 
of mass M* instead of the bare mass M. Equation (0) just carries this information to the 
scalar-baryon coupling fields. 



The Walecka model can be obtained as a particular case of either Lagrangian density by 
making the choice m*(a) = (1 — g a a/M). More generally, recent nonlinear models like those 
presented in Refs. [13], |23|, |24| (which are variations of the DSC model) may be interpreted 
as modified Walecka models where, 

and g* (hereafter we will interpret * as refering to effective coupling constants in the medium) 
is now a function of a related to m*(a) by 

m» = l-s>/M. (8) 

By itself, this establishes a class of models since m*(a) is still a generic real function. In the 
usual Zimanyi-Moszkowski (ZM) model JT31, for example, 



m* ZM {a) = (l + g a a/M) ~\ (9) 



and g* = ga , m* ZM (a). Likewise, the identification of g* and m*(cr) for the other models |£2], £3], 
can be easily done. 

We introduce next a modified version of Model I that we refer to as Model II. We keep 



the generalized factor m*(a) but, following the suggestion given in Ref.|13[, we restrict the 
m*(a) dependence in £ to the fermionic kinetic energy term. To this end we define a modified 
covariant derivative p m *{9uj) = ([m*(a)]^ 1 ip — g^ip) and write the Lagrangian in the form 

C = i> {p m * (&,) — M} ip — -F^F^ U + \mlu^ + X - (d^a - m 2 a a 2 ) . (10) 

A rescaling similar to the previous one, ip — > [m*(a)]^ip, transforms 

{p m *{g u ) — M} — > {p(g*J - m»M} , (11) 

thus generating a new class of nonlinear models. The connection with the Walecka model is an 
extension of the previous one and proceeds through the change, 

CNL = £walecka(g<r ~ > g* a \ 9u ~ > fiC) • (12) 



Here g* is connected with m*(a) by Eq.(|) and 

^ = m*(a). (13) 
9u 

Equations (0) and (|T2|) help to understand the different types of nonlinear couplings showing 
that they are, in fact, effective Walecka models. An important point to emphasize is that given 
the Lagrangian of Eq . (|10|) , the scaling of the effective vector coupling constant in the medium 
is fixed by (0). 

We now explore the possibility of extending the same scaling to the effective scalar coupling 
constant in the medium. Thus, we look for a function m* (a) such that the following constraint 
is satisfied: 

is. = im = m * . (14) 

9a 9w 

It turns out that the only function that fulfills this requirement is m*(a) = m* ZM (a). We call 
the Model II Lagrangian with m*(a) = rn* ZM (a) the modified ZM model (MZM) and is the 
first hadronic model which exhibits this property. It couples the scalar a to the vector u field 
and gives results for the nuclear matter observables substantially different from the usual ZM 
model. A comparison follows. 

We took for the usual ZM model, as given by Eq.flj]), the coupling constants used in 
Rcf. [|l]|. For the MZM model, Eq. ([l2|) , we fit the parameters to the saturation density 
and the experimental binding energy per nucleon to obtain, C% = g^M 2 jm^ = 443.3 and 

= g 2 M 2 /m 2 J = 305.5 |l4j]. In the table below we show, for both models, the numbers 
obtained for the compressibility k, the scalar potential S, the vector potential V (in MeV) and 
the baryonic effective mass m*\ 



Model 


m 


K 


S 


V 


ZM 


0.85 


225 


-141 


82 


MZM 


0.72 


156 


-267 


204 
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Compared to ZM, the MZM model gives, simultaneously, a smaller m* and a smaller k. 
This is unlike what is observed when comparing ZM with the Walecka model. This difference 
in predictions is explained by the nonlinear scalar-vector coupling contained in the MZM La- 
grangian. Furthermore, it is known that the ZM model gives poor results for the spin-orbit 
splitting (which is strongly dependent on the quantity V — S) when used in finite nuclei calcu- 
lations P^ ]. In the MZM model, this quantity more than doubles that of ZM, suggesting that 
MZM may improve upon the former in this particular direction. 

It is also interesting to remark that at low energies the slope of the real optical potential 
(given by 1 — m*) provides information regarding the expected value of m*. The experimental 
value, in the limit of infinite mass number and zero radius, gives m* around 0.6 J23|, |25|. This 
result does not lend support to the ZM model, but tend to favour the Walecka and MZM 
models instead [14]]. Finally, k is smaller in the MZM model than the "empirical" prediction 
k = 210 ±30 MeV. Despite this, the advantage of MZM here is that, unlike the nonlinear a — uj 
model [pq] , it does not present anomalies in the EOS for any value of k. 

We proceed, now, to generalize the ZM and the MZM models. We recall that both originate 
from Eq.(|J) for the particular choice m*(a) = m* ZM (a) . As a consequence of this choice the 
scalar effective coupling constant scales as (g*/g a ) = m*. However, different choices of m*(a) 
can be made meeting the requirement of Eq.(Rj). 



2.1 Model I 

From Eq. (|7]) we generate a family of models by choosing a scaling given by 

Model I: g*Jg a = m* a g*Jg u = l (15) 

The equation of state for such a family of models will be shown in the next section in a general 
expression along with that of Model II. For Model I we consider two particular cases: the 
Walecka model (a = 0) and ZM model (a = 1). 



S 



2.2 Model II 

Defined by Eq.([L2|) and the scaling 

Model II: g*/g tr = m* a gl/g» = m* (16) 

(the MZM model is obtained for a = 1). For a given value of a, the function m*(a) is defined 
by Eq.([U]) together with (|8|) as in the Model I case. 

2.3 The Equations of State for the Models 

When the meson fields in the Lagrangian of both models are substituted by their mean values 
we arrive at the mean field approximation (MFA). In Model I the scalar- meson field a is a 
function of the scalar density (p s ) exclusively In Model II, a depends on p s but also on the 
baryonic density pb- For rotationally and translationally invariant symmetric nuclear matter, 
the MFA equation for the scalar fields reads 



g<j m* a+1 

a 



m 2 M (1 — a)m* + a 



Mp s + p (^L 



2 

* 2 

m p b 



(17) 



with (3 = or 1 for Model I or Model II, respectively. Equation (|17| ) shows clearly how Model 
II and its generalization mix the scalar and vector fields. The scalar and baryonic densities are 
related through, 

Ps (Cl\ ({l-ot)m* + a\ fS\ f V 



77 -/5 



where 



and 



p b ~~ \C*J\ m* 2a+1 - 2 P )\V) \Mm* 

S=- g ;a = -M(l-m*), (19) 



V = (C 2 jM 2 )m^/p b . (20) 

This density ratio weighs the content of each model. Model II has an additional term making 
its contribution to the cx-field, and consequently to the effective mass, larger than that of Model 
I. 



The expressions for the energy density and pressure at a given temperature T can be found, 
as usual, by the MFA average of the energy-momentum tensor. They read, 



and 



where 



C, 2 t2/J 9 M 4 /l-m*\ 2 1 7 /" „ k 2 , 



7 



, r/-'7/(// ; , -//,). (23) 

Here 7 is the degeneracy factor (7 = 4 for nuclear matter; 7 = 2 for neutron matter), nj, and 
stand for the Fermi-Dirac distribution for antibaryons and baryons with exponent (E* ± v) /T, 
respectively. The energy E*{k) is given by 

E*(k) = (k 2 + M* 2 )^, (24) 

while the effective chemical potential which preserves the number of baryons and antibaryons 
in the ensemble is defined by v = p — V , with /1 the thermodynamical chemical potential. 



According the Hugenholtz-van Hove theorem [p7l, the Fermi energy must be equal to the energy 



per baryon at saturation density. Therefore, the following relation has to be satisfied, 

- = V + E*( Po ). (25) 

Po 

We end the section by presenting an analytical general expression for the compressibility 
valid in both models, 

* = + V, + Aw* + 9 {% (wTW^fW + 2 ' 8 WTs) ) ) • (26) 

where kp is the Fermi momentum and all quantities are calculated at nuclear matter saturation 
density p . 
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3 Results and Discussions 



We have conducted calculations with these generalized models at T = 0, for different values of 
a, requiring Ej, = —15.75 MeV at p a = 0.15 fm~ 3 . We begin this section by discussing Model 
I. The question we ask ourselves is whether, by varying a, we can obtain a simultaneous fit for 
k, M* and V — S. The answer to this question is negative. This can be better appreciated by 
looking at Fig. 1 where these quantities are plotted as a function of a. No region in the plot 
gives a simultaneous "reasonable" agreement for the discussed quantities. By "reasonable" we 
mean m* around 0.6 (which, we shall see later, fixes the values of V — S at approximately 680 
MeV) and a compressibility k = 210 ± 30 MeV which is the value of derived from the energy 



of the breathing mode of doubly magic nuclei [p8[| . This conclusion agrees with that of Greiner 
and Reinhard using a generalized "ansatz" applied to the same model. A word about the 



nuclear compressibility; using Fermi-liquid Landau theory, Brown |29| has argued strongly in 
favour of a lower value of k. Since the subject is still under debate we use the value of Ref . [p8f 
as an approximate upper limit. 

The results for Model II are shown in Fig. 2. The compressibility shows a weak dependence 
on a in the region of a > 1, and reaches a minimum close to a = 1 which corresponds to 
the MZM case. Particulary interesting are the results for the region a < 1 where m* and 
the difference V — S show an improvement over the results of the MZM model given in the 
previous section. For example, for a = 0.9 Model II gives M* = 558.3 MeV, V — S = 687.6 
MeV and a compressibility k = 166.2 MeV. It is also interesting to note that for a = 0.88 one 
gets M* = 507.0 MeV, and V — S = 785.1 MeV which are very close to the values obtained in 
the usual linear Walecka model but, this time, with a softer compressibility, k = 181.3 MeV. 

Figure 2 shows also that for a < 0.79 nuclear matter does not saturate. This includes the 
a=0 case, which corresponds to the Walecka model if m*(a) is as in Eq.(|8]). In addition, Figs. 1 
and 2 indicate that M* increases with a thus decreasing the vector potential due to (|2"5|). We 
have allowed for higher values of a in the calculations to reach the curious situation where V 
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vanishes, and model I and II become degenerate. Nuclear matter saturation is achieved with 
just the scalar field. This situation occurs for a ~ 12.8, which is the maximum value that this 
parameter can have (beyond it V becomes negative), and for C% = 315.36. The results for this 
special case are: M* = 885.8 MeV, V - S = -S = 52.4 MeV and K = 65.8 MeV. 

With some insight into the two models, we discuss here if in any M* could give model- 
independent properties of nuclear matter. As we argued before the effective mass is a mani- 
festation of the dynamical relativistic content of any particular model. On the other hand, if 



two models have the same M*, they have the same density ratio as given by Eq.(18). This 
quantity is shown in Fig. 3 as a function of a. We notice that Model II, unlike Model I, can 
acquire a strong relativistic character. Going back to Figs. 1 and 2 we notice that the same 
M* can be obtained in both models for different values of a. Consequently, they will give the 
same relativistic ratio. If the situation is such that two different models give the same M*, the 
question is whether other observables can be determined likewise. For instance, for a given M*, 



the potentials S and V are fixed by Eqs.(|l9l) and (P5|) , respectively. Therefore, the quantity 
V — S and M* are directly correlated and carry the same physical information for each partic- 
ular model. As for the correlation between k and M*, the situation is different having a more 
model-dependent character. 

In order to learn about the compressibility k, we have calculated its individual contributions 
according to Eq. (|26"D . To the first three terms of (|26|) we will refer as k-i and k 3 , respectively. 
An additional fourth term (k 4 ) is present for Model II (/? = 1). These quantities are plotted 
in Figs. 4 and 5. K\ and K2 are completely determined by M* since they depend directly on S 
and V; K3 and K4, however, show an extra dependence on the slope of M*(p) which is negative 
and carries information about the scalar field. For Model I, Fig. 1 shows that as M* increases 
the compressibility k decreases. This explains why, in the Walecka model where M* is small, 
k is so high. In the usual ZM model the opposite is true. In Model II the changes in K3 and K4 
are such that despite the increase in M* the total compressibility remains almost constant. 

Finally, in Fig. 6 we show the effective mass m* — solution of Eq.(Rf) — for different values 
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of a. 



4 Conclusions 

In summary, we have shown the equivalence between two descriptions of derivative scalar cou- 
pling models which we now interpret as effective Walecka models. In them the effective coupling 
constants depend on the density and are completely determined by the effective nucleon mass 
m*. The first model has only one effective coupling g* related to m* through m* = 1 — g*a/M. 
The second includes also an effective vector coupling g* which is always given by g* = m*g LJ . 
We have also shown that for a particular choice of m* = m* ZM (a) both effective meson coupling 
constants scale with m*. We have generalized this scaling by introducing an extra parameter 
a and concluded that, in a theory with three free parameters, Model I can not provide good 
results for the level splitting in finite nuclei and at the same time adjust the experimental k 
and m*. We showed, however, that it is possible to choose the parameter a so as to give values 
of m* , k and V — S similar to those obtained from the nonlinear a — u model. We conclude 
that changing only the scalar coupling (Model I) definitely improves k but with no hope of re- 
producing the spin-orbit splitting for finite nuclei, thus supporting the conclusions of Ref.|22 . 



This is a direct consequence of the large values M* that Model I gives when k goes in the right 
direction as a changes. Model II, instead, through the inclusion of a mixed coupling between 
the scalar and the vector fields may be the way to expect improvement in the calculation of 
the observables. 
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Figure 1: The nucleoli effective mass (M*), the difference between the vector and scalar 
potentials (V — S) and the compressibility k as function of a for Model I. 



Figure 2: The nucleon effective mass (M*), the difference between the vector and scalar 
potentials (V — S) and the compressibility k as function of a for Model II. 



Figure 3: The relativistic ratio between scalar and baryonic densities (p s / p&) as a function of 
a for both models. 



Figure 4: The components k±, k 2 and Ks of the compressibility k as a fuction of a for Model 
I. 



Figure 5: The components k±, k,2, ^3 and k 4 of the compressibility k as a fuction of a for 
Model II. 



Figure 6: The effective nucleon mass as a function of the scalar sigma field (u = g a a/M) for 
different values of a. 
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